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Introduction 


Justification du choix et Eclaircissement de la 
problématique 


Les réseaux deau potable sont de plus en plus souvent équipés de sys- 
tèmes de télégestion des installations et des commandes. Laccroisse- 
ment de ces équipements accompagne généralement une augmenta- 
tion de la complexité du système de distribution. 

En effet, pour des raisons de sécurité mais aussi de rentabilité, beau- 
coup de réseaux résultent de l’interconnexion de sous-réseaux qui 
étaient auparavant gérés de façon autonome. Ceci augmente la diffi- 
culté de choisir, au cours d'une journée, les consignes à donner aux 
différents organes d'un réseau deau potable. La complexité de la ges- 
tion vient aussi de la prise en compte par les gestionnaires des réseaux 
de deux objectifs contradictoires : la sécurité d'approvisionnement et 
la minimisation des coûts de production et de transport de l'eau. Pour 
aider les opérateurs dans la gestion des réseaux, de nombreux out- 
ils d'automatisation ont été conçus depuis une vingtaine d'années. Ils 
peuvent être regroupés en deux catégories : ceux faisant appel à des 
systèmes experts s'inspirant initialement des méthodes utilisées par 


les hydrauliciens, et ceux reposant sur des méthodes doptimisation. 


Abstract 


This paper presents a 24 hours optimization models of a drinking 
water network. The method used is based on linear programming. 
This method allows obtaining optimal commands in real time. In 
this paper, we develop two points : the demand forecast on an area 
and the taken into account stochastic aspect of that demand in the 
optimization models. 


Valeur ajoutée et solution proposée 


Notre objectif étant de proposer des méthodes doptimisa- 
tion pour des réseaux de taille moyenne, nous avons choisi les 
méthodes doptimisation demandant peu du travail dexpertise 
préalable sur le réseau étudié. Dans la première partie, nous 
identifierons les principales variables et nous présenterons le 
modèle général doptimisation dit « déterministe ». Il sagit 
d'un modèle où la consommation est considérée comme conue 
Nous présenterons ensuite les modèles de prévision de la de- 
mande permetant d'affiner le paramètre de consommation dans 
les modèles doptimisation. Enfin, nous verrons comment nous 
pouvons modifier les modèles doptimisation pour tenir compte 
de l'incertitude de prévision. 


Explication de la méthode et 
de l'algorihme 


Le problème d'optimisation 


Le but des modèles que nous construisons est dobtenir les commandes à passer sur un réseau 
deau potable afin de minimiser les coûts de prodution et de transport de l'eau. Les organes de 
commande peuvent être des pompes, des vannes ou encore des stations de traitement de l'eau. 


Les variables utilisées 


Soit k le nombre dorganes présents sur le réseau que nous optimisons.On appelle commande 
élémentaire, un vecteur à k composantes traduisant l'ensemble des commandes du réseau. Il 
y a alors au plus 2^k vecteurs de commande élémentaire comportant chacun k composantes 
valant 0 ou 1. Nous précisons qu il y en a au plus 2k car certains vecteurs de commande po- 
tentiels ne sont pas réalisables d'un point de vue physique. Cest le cas, par exemple, de toute 
commande consistant à mettre en route une pompe débitant dans des canalisations où toutes 
les vannes seraient fermées. Les différents vecteurs de commande seront notés par la suite 

. À titre d'exemple nous pouvons écrire les vecteurs de commande élémentaire d’un réseau, 
représenté figure 1, et comportant une station de traitement et deux vannes. 


v Ss adh L= Exemple d'un réseau comportant trois organes de com- 
mande. 
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Nous pouvons alors détailler le cas du vecteur K6 qui traduit le cas suivant : 

l'usine de traitement deau fonctionne ; 

ela vanne d'alimentation de R1 est ouverte ; 

ela vanne d'alimentation de R2 est fermée. 

En fonction de la configuration hydraulique de ce réseau, certains des vecteurs seront à écarter. Nous in- 
troduisons des variables ai(t) qui représentent la fraction de temps d'application de chacun des vecteurs de 
commande pour chaque pas de temps [t- 1;t] De cette définition découle la propriété : 0 < ai(t) < 1. On peut 
alors exprimer la commande globale sur le réseau par : 


K (t)= $a (t) K 
Le modèle déterministe 


À partir de ces définitions, nous pouvons bâtir un modèle d'optimisation en programmation 
linéaire. Il consiste en la minimisation d'une fonction coût sous plusieurs contraintes phy- 
siques. Nous notons Qri(t), les débits représentant le volume deau entré dans le réservoir r au 
cours du pas de temps|t- 1;t] si la commande i est appliquée et ni(t) représentant les coûts de 
cette commande i. De plus, nous noterons Cr(t) la consommation de la zone dépendant du 
réservoir r au cours du pas de temps[t- 1;t].La mise en œuvre du modèle suppose de calculer, 
avant loptimisation, les débits et les coûts. Pour le calcul des débits, nous pouvons soit utiliser 
les données enregistrées aux entrées et sorties des réservoirs sur le réseau soit une modélisa- 
tion à l'aide du logiciel Zomayet développé par l'unité « Ouvrages et réseaux hydrauliques » 
du Cemagref. Nous présentons le modèle dans le cas général où il y a T pas de temps. Dans la 
pratique nous travaillerons au pas de temps de l'heure et sur un horizon de 24 heures, nous 
aurons donc T = 24. Le modèle sécrit de la façon suivante : 


tmn = $a, (6), (0) 

(Y8 e (1;T -1},Yre (1;R},V, 0È $a (t)Q,; O-Èc, (E)E Vrm * 
(2)V0 e€ (1;T —1},Yr e (1;R},V, CDD ()9, O-Ÿc, ()2 Vin 
OVre(HR).V, (0)+ È $a (90, ()- ÈC, (= Von 


(48E (T), Ya (8)=1 


VOe(IT-1),Vie(1;1),a,(8)>0 


Les contraintes de type (1) et (2) traduisent les limites maximales et minimales dans chacun 
des réservoirs au cours de l'horizon doptimisation. Elles sexpriment grâce à un bilan entrée/ 
sortie dans chaque réservoir et à chaque pas de temps. Les paramètres Vrmax et Vrmin tra- 
duisent les volumes utilisables lors de l'optimisation La différence (Vrmax- Vrmin) peut se 
calculer comme le volume physique du réservoir auquel on retranche le cas échéant la réserve 
incendie et les marges de sécurité choisies par le gestionnaire du réseau. La contrainte (3) 
permet de fixer un niveau final qui consiste à remplir les réservoirs. Enfin, la contrainte (4) 
découle de la façon de construire les variables ai(t). 


La prévision de la demande 


Nous allons présenter ici la prévision de la demande des zones de consommation, dans le 

but d'améliorer les performances du modèle doptimisation. Cette prévision se fera au pas de 
temps de l'heure alors que beaucoup de modèles de prévision de la demande deau potable ou 
délectricité (Canu, 1990 ; Mangeas, 1993 ; Jain, 1994) travaillent sur un pas de temps journal- 
ier. Pour construire les modèles de prévision de la demande nous avons utilisé les données 

du syndicat intercommunal des eaux de la région de Sarralbe qui se situe à la limite entre les 
départements de la Moselle et du Bas-Rhin. Le réseau de ce syndicat est affermé à la Com- 
pagnie générale des eaux. Nous avons comparé deux grands types de modèles : les réseaux de 
neurones et les séries temporelles. Nous ne présentons ici que les modèles construits à l'aide 
des réseaux de neurones, la différence des résultats obtenus entre les deux outils nétant que 
peu significative. 

Lobjectif initial était de bâtir des modèles prédisant la demande d'une zone pour les 24 heures 
à venir conformément aux besoins du modèle doptimisation présenté ci-avant. Pour atteindre 
cet objectif, nous avons tout d'abord construit un modèle déterminant la demande de l'heure 
à venir. Nous avons ensuite établi un modèle pour la prévision de la demande sur plusieurs 
heures consécutives. 


Les données utilisées : application au syndicat 
intercommunal des eaux de la région de Sarralb 


Nous travaillons sur une partie du réseau de ce syndicat intercommunal pour lequel nous 
disposons de données heure par heure allant de juillet 1993 à décembre 1997. Ces différentes 
données sont issues de compteurs situés à l'entrée et à la sortie des réservoirs ainsi que des 
enregistrements des hauteurs deau dans ces mêmes réservoirs. Nous avons ainsi pu isoler une 
zone de demande, appelée « zone ouest », présentée figure 2, représentant environ 10 000 hab- 
itants pour une consommation journalière moyenne de 700 m3 

avec des journées pouvant atteindre 1 100 m3 

en pointe. Des zones de cette taille comportent une variabilité horaire et journalière impor- 
tante, ce qui rend la prévision plus difficile quen zones urbaines. 
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À Figure 2. - Schéma du réseau du SIER de Sarralbe. 


Dans un premier temps, nous avons éliminé, à l'aide de règles de validité, les périodes com- 
portant des fuites importantes ou celles mêlant de la consommation à du transit entre rés- 
ervoirs. Une fois ces opérations de tri effectuées, nous pouvons travailler sur 1 025 journées 
complètes. 


La construction des modèles 


Le choix des entrées 

Pour la construction de modèles de prévision de la demande, plusieurs auteurs (Canu, 1990) 
ont utilisé des données météorologiques ou climatiques. Nous avons choisi de nous affranchir 
de ce type de données et ceci pour quatre raisons : 

1. Les liens entre des données climatiques et la consommation semblent difficiles à établir car 
il faut choisir le facteur explicatif adéquat (nombre de jours consécutifs de sécheresse, nébu- 
losité faible, seuil de température..). 

2. Lutilisation de prévisions météorologiques pose le problème des incertitudes du modèle 
météorologique qu'il faudrait transposer dans le modèle de prévision de la demande, ce qui 
est délicat à maîtriser. 

3. Les modèles que nous construisons devront pouvoir être implantés sur un site central où le 
réseau est géré en temps réel. Les données dentrée du modèle devront être facilement accessi- 
bles ce qui nest pas le cas pour des données météorologiques. 

4. Les données météorologiques sont implicitement contenues dans les consommations 
passées. Par exemple, lorsque la prévision de demande se fait au cours d'un épisode clima- 
tique chaud et sec, les consommations des heures et des jours précédents traduiront cet ép- 
isode par des valeurs élevées. Nous choisirons d'utiliser plutôt ces consommations passées 
dites données endogènes que les données climatiques. 

Les variables exogènes que nous avons alors recensées pour figurer en entrée du réseau de 
neurones et qui sont apparues comme significatives sont le mois de l'année, le jour de la se- 
maine, l'heure de la journée et le fait que le jour considéré, corresponde à une période de 
vacances scolaires (oui ou non). Nous aurons besoin dans l'optimisation d'un réseau deau 
potable de connaître la demande sur plusieurs heures consécutives. Idéalement les modèles à 
bâtir devraient prévoir la demande sur 24 heures ce qui correspond à l'horizon doptimisation 
que nous utilisons. Dans les faits, notre optimisation débutant à 6 heures du matin avec len- 
semble des réservoirs pleins, nous savons qu'il nous faut une prévision fiable de la demande 
un peu plus tard dans la journée. Cest pourquoi nous avons examiné un modèle prévoyant 
les demandes sur 12 pas de temps de 13 heures à 24 heures. Le modèle utilisé est représenté 
figure 3. 

En plus du biais, ce réseau comporte 30 entrées réparties de la façon suivante : 

6 neurones, nommés de DIMANCHE à VENDREDI, prenant la valeur 0 ou 1 pour 

définir le jour de la semaine (si tous ces neurones valent 0, il s'agit d'un samedi) ; 

11 neurones définissant de la même façon le mois de l'année ; 

12 neurones appelés Di qui prennent la va leur de la demande des 12 premières heures de la 
journée en m3 ; 

e] neurone valant 0 ou 1 pour préciser si lon se place lors d'une journée de vacances scolaires. 


De plus ce réseau de neurones comporte une seule couche cachée avec 2 neurones (en 
plus du biais) au sein de celle-ci. Enfin, nous complétons le réseau par une couche de 
sortie composée de 12 neurones qui vont nous permettre dobtenir les 12 demandes de la 
deuxième partie de la journée. 


Les calculs 


Nous représentons figure 4 le coefficient de corrélation entre les valeurs attendues et les 
valeurs prédites pour un calcul débutant à 12 heures et ce sur la partie des données ayant ser- 
vi à la validation du modèle. 
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À Figure 4. - Coefficients de corrélation du 


À Figure 3. - Réseau de neurones pour la modèle sur plusieurs heures. 


prévision sur un horizon de plusieurs heures. 


La prise en compte de l'aspect 
stochastique de la demande 


Plusieurs auteurs se sont intéressés à la prise en compte de l'aspect stochastique de la de- 
mande dans l'optimisation des réseaux deau potable. On peut en particulier citer les travaux 
d'Aleksandrov (1984), qui s'intéressent également à un système de distribution deau potable, 
mais avec des choix déchelles spatiales et temporelles très différentes des nôtres. Pour notre 
part nous sommes partis des travaux présentés par Golstein (1973), qui fut l'un des premiers 
à introduire l'incertitude de la demande dans un modèle de distribution tout en conservant 
des modèles de type « programmation linéaire ». Nous proposons, en fait, deux méthodes 
pour prendre en compte le caractère aléatoire de la demande dans les modèles doptimisation. 
Avant de présenter ces méthodes nous allons définir les paramètres utiles dans la prévision de 
la demande 

Les modèles de prévision 

Les modèles de prévision de la demande que nous avons obtenus peuvent 

sécrire : 

Cr(t) = fr(P(t), y(t)) + er(t) 

Le terme fregroupe toutes les variables endogènes et le terme ytoutes les vari- 

ables exogènes utilisées pour la construction du modèle de prévision. Le terme 

er(t) représente l'erreur de prévision. Cest pour tenir compte de ce terme que- 

nous voulons introduire la programmation stochastique. La distribution des 

termes er(t) est aléatoire de moyenne nulle et décart-type donné. 


On peut, bien entendu, regrouper les deux termes et présenter alors la demande Cr(t)comme 
ayant une distribution gaussienne de moyenne u (Cr (t)) et décart-type o (er). 

La distribution de la demande 

Pour chaque variable aléatoire Cr(t), nous noterons u(Cr(t)) la moyenne et o(Cr(t)) = o(er) 
lécarttype. Si nous supposons de plus que ces demandes sont statistiquement indépendantes, 
les sommes partielles.c, = ase (t) 


t=1 


suivent également une | ar gaussienne de moyenne H(C,e )= D u(c,(r)) 


et decarttype o(c)=/Xo"(c,(r)) 


Nous aurons besoin de discrétiser la distribution de la demande. Dans ce cas, nous utiliserons 
les propriétés d'une distribution gaussienne rappelées figure 5. 
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À Figure 5. - Forme générale d'une distribution gaussienne. 


prévision de la demande. La première dite « formulation serrée » nous permet denvisager 
les cas les plus extrêmes de demande. Ces cas extrêmes sont soit des demandes très faibles 
entraînant un risque de débordement des réservoirs soit des demandes très fortes entraînant 
un manque deau dans certains réservoirs. La seconde méthode proposée sera celle des 
«problèmes àcontraintes probabilisées » qui permet dese prémunir contre un risque fixé par 


l'utilisateur du modèle. 
Formulation serrée 


Cette formulation correspond au cas où la violation d'une des contraintes aurait des 
conséquences très fâcheuses. Le programme est alors constitué par des vecteurs qui satisfont 
les contraintes dans tous les cas, cest-à-dire pour toute probabilité q parmi les états de la na- 
ture possible Q 


Les contraintes du problème peuvent alors 
sécrire : pour tout q appartient a Q, A(q).X < B(q), X > 0 
Si l'ensemble des couples [A, B] est fini, la résolution du problème stochastique se ramène au 


problème déterministe : 
ax E(C,X ) sous 


Vae(lk),A(q).X <B(q) 

X 20 
Le domaine des admissibles est alors un polyèdre éventuellement non-vide limité par les k 
contraintes précédentes. Dans la pratique, nous pourrons considérer que l'ensemble des cou- 
ples [A, B] est fini en perdant l'information sur les queues de distribution. 
Nous noterons : 
e SrO le nombre de pas de discrétisation pour chaque réservoir et à chaque pas de temps ; 
e pour s appartient a [1; SrO], Cr&s, la valeur discrète de la demande 
aléatoire cumulée. 
Dans la pratique nous choisissons Sr0= 5 avec des valeurs extrêmes rendant négligeables les 
queues de distribution. Nous pourrons choisir, par exemple, la moyenne + ou - 4 fois l'écart- 


type. 


Écriture du programme 


Znin = > D (0), (8 
geal j=l 


(Vee(LT-1},Vre(LR},Vse(ls,,},V, (0)+ D Şa (t) 2n (t)-V, ma E Cros 


i=l i=l 


(2)V8e(LT-1), Vre (1; R}, Vse (1; s), V, (0)+ Sa (t)Qu (E )-V an > 


=] i=] 


(3)vre (1; R}, Vse {Ls} V, LHÈŸe (1)Q,()-Vn = Cr 


(4)V8e (ET), Ÿ @ (8)=1 
vVae(lT-1},Vie(l:1),a(8)z0 


Contraintes incompatibles 


Il est possible que cette formulation aboutisse à un problème non réalisable. Le domaine des 
admissibles constitué par les différentes contraintes peut, en effet, être vide. Dans notre cas, 
les contraintes de type (3) sont incompatibles entre elles. Ceci est dû à la présence dégalités 
Nous aurons donc à transformer ces contraintes. 


Les problèmes à contraintes 
probabilistes 


r, , 
Les contraintes du problème sécrivent: Wie (1, m)  P Xa, x, < b, |2 p,,0< p, Sl. 
j=l 


Ceci signifie que le programme vérifie la ième inégalité avec une probabilité P au moins égale 
à pi. Pour la résolution de ce type de problème on utilise lľéquivalence suivante : 


Vie (1,m),P Şax, <b, |> p,,0< p, S1. 


j=l 


Dans laquelis b est déterminé par: Í g. (b )db, =p, et Q, (b ) 


b, 


1 m 1 T 


est la densité de distribution de la composante bi. Dans le cadre d'une distribution gaussienne 
la densité de distribution est bien connue et il est aisé d'inverser la relation précédente. Par 
exemple, si pr8= 0.9, cest-à-dire si la probabilité de débordement du réservoir r au pas de 


temps [ 0- 1; OÜJest de 90 %, alors: Ca = u (C, -1.41x0 (C, ))- 


Application et résultats 


Nous allons présenter les deux méthodes sur un exemple. Nous pouvons les mettre en œuvre 
grâce au logiciel doptimisation Cplex de la société Ilog. 


Le réseau exemple 

Nous utilisons l'exemple présenté figure 1. Formulation serrée Le temps de calcul est de 0,5 
sec CPU. La valeur de la fonction objectif servira de majorant pour la méthode suivante. Nous 
noterons 100 cette valeur. 

Formulation avec des contraintes probabilistes 

Pour s'assurer de la cohérence de nos calculs, nous pouvons construire un modèle pour lequel 
toutes les probabilités valent 99,997 %. Il nous faut là aussi relâcher les contraintes de stock 
maximal (contraintes de type (1)) dans les derniers pas de temps de l'horizon doptimisation. 
Ceci donne la même valeur de la fonction coût que dans le cas de la formulation serrée. 

Une fois cette vérification effectuée, nous modifions les valeurs des différentes probabilités 
Nous présentons le résultat des différentes simulations dans le tableau suivant. 


Pe = Po = 99,997 % 


p? = 99,997 % | Contrainte (1) utilisable Contrainte (1) utilisable 
de 1 à 20h de 1à22h 


Contrainte (1) utilisable Contrainte (1) utilisable 
de 1 à 22h de 1 à 22 h 








Dans ce tableau, nous comparons les coûts des solutions pour lesquelles les probabilités des 
contraintes de type (1) et (2) sont égales pour tous les pas de temps et valent soit 90 % soit 
99,997 % ainsi que les coûts des solutions pour les mêmes valeurs de probabilité pour les 
contraintes de type (3). Nous notons également, pour chaque case du tableau, l'intervalle 

de temps sur lequel les contraintes de type (1) sont utilisables. Ce dernier résultat résulte du 
problème d'incompatibilité présenté ci-avant (voir Contraintes incompatibles). 

Nous pouvons en conclure que le choix de la probabilité de remplissage des réservoirs (con- 
trainte de type 3) a une influence forte sur le coût de la solution. Le choix des probabilités 
liées aux niveaux minimum (de type 2) et maximum (de type 1) dans les réservoirs ne modi- 
fie pas les coûts dans notre exemple mais peut éventuellement permettre de rencontrer moins 
de problèmes d’incompatibilité. Les temps de calcul de ces différentes simulations sont égale- 
ment de l'ordre de 0,5 sec CPU. 


Des perspectives pour un élargissement 


Les deux méthodes proposées permettent de tenir compte de l'aspect aléatoire de la demande 
en eau potable après la mise en place doutils de prévision de cette demande. Dans les deux 
cas, le nouvel outil de gestion en temps réel que nous proposons est plus sécuritaire car il tient 
compte de la distribution de la demande au lieu de n'utiliser que la valeur moyenne comme 
nous le faisions jusquà maintenant. 

Dans la formulation serrée, nous obtenons une solution très sécuritaire mais qui entraîne des 
surcoûts importants. Dans la formulation probabilisée, le gestionnaire du réseau peut choisir 
les différentes probabilités en fonction des réservoirs et des moments de la journée. Cest alors 
lui qui maîtrise la sécurité de l'outil de gestion en temps réel. Il devient ainsi possible d'affiner 
les marges de sécurité forfaitaires que les gestionnaires de réseaux utilisent. 

Au cours de l'étude (Guhl, 1999), nous avons intégré plusieurs contraintes physiques nouvelles 
dont la résolution na pas été présentée dans cet article. Il s'agit d’une part de ce que nous 
appelons les « contraintes de continuité ». Elles sont liées à l'obligation, sur certains ouvrages 
comme des stations de traitement de l'eau, d'avoir un fonctionnement par plages continues. 
D'autre part, nous avons construit des méthodes permettant de tenir compte de l'influence de 
état du réseau sur les transits entre ouvrages sans utiliser la programmation dynamique ou 
non-linéaire. 


Ceci nous permet d'appliquer les modèles doptimisation à des réseaux pour lesquels la var- 
iation de hauteur d'eau dans les réservoirs est plus importante que la variation des altitudes 
comme cest le cas dans des régions où la topographie est peu variable. 

La phase de validation des méthodes proposées sur des réseaux deau potable télégérés est en 
cours actuellement. Dans un premier temps nous passons par une phase doptimisation que 
nous qualifions d’« hors ligne » qui permet de proposer aux gestionnaires d'un réseau des 
stratégies optimales de commandes d'un réseau sans les appliquer directement. Les stratégies 
proposées peuvent, en particulier, prendre la forme de courbes consigne dans les réservoirs qui 
remplacent les seuils d'arrêt/démarrage classiques par des valeurs modulables sur 24 heures. 


Application au Maroc 





Le secteur de l'eau potable et de l'assainissement au Maroc a vu depuis les 

dernières décennies, une augmentation très importante de l'accès à leau. 
į , ° $ ° 7 3 7 e, 7 

ce qui rend loptimisation du réseau deau potable au Maroc une nécessité 

voire une urgence. 

l'étude faite dans ce projet peut être considérée comme une opportunité 

pour le Maroc pour amélio Wnnaît des temps 


Bureau d'etude vs 
une grande société 


Si ce projet senvisagerait dêtre réaliser au Maroc, il doit être projeté sur 
7 ° ° e, 7 2 ° b} 
une échelle nationale pour sentir son efficacité. alors cest un travail d'une 

grande entreprise quiļ'est mise en jeu.E 
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